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Abstract 

In this paper a theoretical model of a magnetic white dwarf is stud- 
ied. All numerical calculations are performed under the assumption of 
a spherically symmetric star. The obtained equation of state is stiffer 
for the increasing value of magnetic field \ (B\). The numerical values 
of the maximum mass and radius are presented. Sensibility of the re- 
sults to the choice of the magnetic field strength is evident. Finally the 
departure from the condition of isothermality of degenerated electron 
gas in gravitational field is discussed. 



1 Introduction 

The possibility of the existence of white dwarfs with large magnetic fields 
was predicted by Blackett in 1947. Ever since Kemp et al. in their paper 
presented results of the first detection of magnetic white dwarfs, these kind of 
stars have been intensively studied (Kemp et al. 1970; Putney 1995; Reimers 
et al. 1996; Jordan 1992; Chanmugan 1992; Ostriker et al. 1968; Mathews 
et al. 2000). However, only a small minority (the exact percentage remains 
uncertain) of white dwarfs have detectable magnetic fields. The study of 
magnetic properties of white dwarfs are expected to provide information on 
the role played by magnetic fields in the stellar formation and the main 
sequence evolution. Magnetic fields in white dwarfs are interpreted as the 
relic of magnetic fields of their progenitors. Thus the evolution of the main 
sequence stars with fields of only a few Gauss leads to white dwarfs with 
weak magnetic fields which remain below the threshold detection. Magnetic 
Ap and Bp stars are considered themselves as progenitors of magnetic white 
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dwarfs with the field strength > 100 MG. 

The aim of this paper is to show how the assumption of finite tempera- 
ture and nonzero magnetic field influences white dwarf parameters. In the 
presence of magnetic field the magnetic pressure which corresponds to the 
Lorence force has to be taken into account. However, its inclusion causes the 
deformation of a star from spherical configuration, (Konno et al. 1999). In 
this paper all calculations are performed with the assumption that the aver- 
aged value of the energy momentum tensor gives us the condition for small 
deformation of a star and allows to treat a star as a spherically symmet- 
ric one. This approximation is connected with the existence of the critical 
density pcrit below which a star became deformed. The gained results allow 
to construct the mass-radius relation for a white dwarf with the magnetic 
field and at finite temperature which has been chosen to be equal 10^ K. 
The description of the interior of a white dwarf has been made in terms of 
degenerate, uniform, isothermal, electron gas. Using the equation of state 
of classical gas the atmosphere model has been formulated. The atmosphere 
extent, luminosity, effective temperature and the value of the critical den- 
sity pcrit as functions of the temperature, density and magnetic field strength 
have been examined. The determination of the critical point allows to neglect 
the deformation of the star and perform all calculations in the spherically 
symmetric approximation. 
This paper is organized as follows. 

Section 2 outlines the considered model and contains assumptions which are 
indispensable in determining the form of the energy momentum tensor and 
the value of the critical density Pcrit- In section 3 the numerical results and 
the discussion are presented. 

2 The theoretical model of white dwarfs in the mag- 
netic fields 

This paper presents the theoretical model of a white dwarf in which the 
main contribution to the pressure which supports the star against collapse 
comes from ultra-relativistic electrons . The electrons are supposed to be 
degenerate with arbitrary degree of relativity a; = ^ (in this paper the 
units in which c = h = 1 were used ). The value of x is a function of 
a density. Ions provide the mass of the star but their contribution to the 
pressure is negligible. The complete and more realistic description of a white 
dwarf requires taking into consideration not only the interior region of the 
star but also the atmosphere. The isothermal interior of a white dwarf is 
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completely degenerated and covered by non- degenerate surface layer which is 
in radiative equilibrium. Thus the white dwarf matter consists of electrically 
neutral plasma which comprises positively charged ions and electrons. 
The Lagrangian density function in this model can be represented as the 
sum 

= >Ce + >Cg + >Cqed, (1) 

where C^, Cq describe the electron and gravitational terms, respectively. 
The >Cqed is the Lagrangian density function of the QED theory 



The electron part of the Lagrangian is given by 

where is the covariant derivative defined as = — ieA^ and e is 
the electron charge. The Dirac equation for electrons obtained from the 
Lagrangian function C has the form 

{iji'Di, - m,)il^e = 0. 
The vector potential is defined as An = {Aq = 0, ^i} where 

Ai = ——eHraX^B^. 

The symmetry in which uniform magnetic field B lies along the z-axis has 
been chosen 

The energy momentum tensor is the sum of the perfect fluid part TJ^^** 

r;^^" = {£ + P)u^^uu + 9^^uP, 
and the electromagnetic one 

and can be calculated taking the quantum statistical average 

J-IXV — ^ J-flV -1-^11/ "T J- ^jj 
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In the cartesian coordinates the electromagnetic part of the energy momen- 
tum tensor in the flat space-time has the form 
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and presents the example of the anisotropic pressure. In general the pressure 
of the electron gas in the magnetic field can be written as the sum of two 
parts 



P = Pe + i'QED 



where Pe 
.1 



is the isotropic part of the pressure, Pqed can reach the vahie 
mses that the pressure became anisotn 
for polar coordinates can be written as follows 



±iP^ and causes that the pressure became anisotropic. On the other hand 
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Averaging over all directions allows to obtain the isotropic form of the energy 
momentum tensor relevant for the spherically symmetric configuration 
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The assumption of the positivity of the total pressure P has been made 

1 9 
Pe - -P > 0. 



Any negative contribution to the pressure reduces it, gravity can not be 
compensated and this diminishes the radius of the star. The described above 
averaging approximates a star as a spherically symmetric object with the 
radius R which is smaller than the equatorial radius of the anisotropic star. 
The case of an anisotropic star is calculated in details in the paper by Konno 
et al. (Konno et al. 1999). The condition of the positivity of the total 
pressure establishes a limiting value of the density noted pcrit or the critical 
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point Xcrit and puts limit on the spherically symmetric approximation (the 
deformation of the star is neglected). For densities smaller than the critical 
one the total pressure became negative and there is no gravity compensation 
in this direction. 

The equations describing masses and radii of white dwarfs are determined 
by the proper form of the equation of state. The aim of this paper is to 
calculate the equation of state for the theoretical model of magnetic white 
dwarfs with the assumption of finite temperature. 

Properties of an electron in the external magnetic field have been studied, 
for example by Landau and Lifshitz in 1938 (Landau et al. 1938). In this 
paper the effects of the magnetic field on an equation of state of a relativistic, 
degenerate electron gas is considered. The motion of free electrons in the 
homogeneous magnetic field of the strength B in the direction perpendicular 
to the field is confined by the oscillatory force determined by the field B and 
is quantized into Landau levels with the energy nB, n = 0, 1, . . .. In the case 
of non-relativistic electrons their energy spectrum is given by the relation 

where n = j + ^ + Sz, the cyclotron energy lOc = eB/nie, Pz is the momentum 
along the magnetic field and can be treated as continuous. For extremely 
high magnetic field the cyclotron energy is comparable with the electron 
rest mass energy and this is the case when electrons became relativistic. 
Introducing the idea of the critical magnetic field strength Be = ml/ \ e | 
which equals 4.414 x 10^'^ G, it is easy to distinguish between non-relativistic 
and relativistic cases, for B > Be the relativistic dispersion relation for the 
electron in the magnetic field is given by 



£;n,p, = VpI + + 2eBzn (3) 

where n is the Landau level, pz is the electron momentum along z-axis and 
rUe is the rest mass of the electron. Along the field direction a particle motion 
is free and quasi-one-dimensional with the modified density of states given 
by a sum 

^'"^ 2y\2-6,o] f^dpz, (4) 



6no denotes the Kronecker delta (Mathews et al. 2000). The spin degen- 
eracy equals 1 for the ground (n = 0) Landau level and 2 for n > 1. The 
redefined density of states makes the distinctive difference between the mag- 
netic and non-magnetic cases. The equation (|3|) implies that for n = 0, 
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-^0 = + m| whereas for n > 1 En = \/vi + . These relations indi- 
cate that the quantity m* defined as m*^ = mg + 2eB^n can be interpreted 
as the effective electron mass which is different from the electron mass me 
for n > 1. The number density of electrons at zero temperature is given by 
the following relation 

F 



[2 - (5no]2eB,p^ 



n=0 

The maximum Landau level is estimated from the condition )^ > 0. One 
can define the critical magnetic density /9b 

/9B = 0.802 y-^-f-l'^g cwT^ (5) 

(7 = B^jBf. means a dimensionless parameter which expresses the magnetic 
field strength through the value of the critical field for densities lower 
than Pb only the ground Landau level is present. 

Assuming the finite temperature of the system which affects the electron 
motion in the external magnetic field the number density of electrons can be 
expressed now as 

^e = ^E:r=o[2-Mx (6) 

(/o,o,+(^A, 1 + 27n) - /o,o,-(^A, 1 + 27n)) 
where the Fermi integral 

h,r,^{u,a) = J (7) 

has been used. The fj,^, and denote the chemical potential and temperature 
in the star interior and z = n-i^/rrie, t = ^BT^/me, u = z/t (Lai 2000; 
Bisnovatyi-Kagan 2000). Finite temperature together with decreasing value 
of the strength of the magnetic field tend to smear out Landau levels. The 
introduction of the concept of Fermi temperature 

Tf = Ep/kB = (me/A;B)eF (forp < ps), (8) 

and the magnetic temperature 

Tb = ^ (9) 

= ^ (VI + 2rama^7 + 27 - VI + 2nmax7) 
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where AEb is the energy difference between the riL = nmax level and the 
iT'L = "-max + 1 level, is the convenient way for describing how the properties 
of free electron gas change at finite temperature and at the presence of the 
magnetic field. For p = pb these temperatures equal Tp = Tb- The influ- 
ence of the magnetic field is most significant for p < Pb and T < Tb when 
electrons occupy the ground Landau level. In this case one can deal with 
the strong quantizied gas and the magnetic field modifies all parameters of 
the electron gas. For example, for degenerate, non-relativistic electrons the 
pressure is proportional to p^ and this form of the equation of state which 
one can compare with that for the case B=0 for which P ~ p^^^. When 
p > Pb the Fermi temperature is still greater than the magnetic tempera- 
ture Tb. Electrons are degenerate and there are many Landau levels, now 
the level spacing exceeds /cbT. The properties of the electron gas are only 
slightly affected by the magnetic field. With increasing temperature, there is 
the thermal broadening of Landau levels, when T > Tb, the free field results 
are recovered. For T 2> Tb there are many Landau levels and the thermal 
widths of the Landau levels are higher than the level spacing. The magnetic 
field does not affects the thermodynamic properties of the gas (Lai 2000 ). 
The total pressure of the system can be described as the sum of the pres- 
sure coming from electrons and ions plus small corrections coming from the 
magnetic field 

P = Pe + PQED. 

The contributions from the same constituents form the energy density 

e = Cion + EQED, 

where Pqed = eqed = ^B^- The electron pressure is defined with the 
use of the Fermi integral 

n = ^Er=o[2-Mx (10) 
/_i,2,+(zA, 1 + 2jn) + 1-1,2 A^/t, 1 + 2jn)). 

The ion energy density has the simple form 

£ion = neTTiB, (H) 

where uib is the baryon mass. The ions are very heavy and are treated as 
classical gas. The influence of the magnetic field for ions is very small and 
this correction to the equation of state has been neglected. 
The pressure and energy density dependence on the chemical potential 
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determines the form of the equation of state, which is calculated in the flat 
Minkowski space time. The obtained form of the equation of state is the 
fundamental input in determining macroscopic properties of a star. Outside 
the static spherically symmetric body the metric has the form 



e-i^)dt^ - e^Wdr^ - r^dO^ - sin^ Odip^, 



where z/(r) and A(r) are metric functions. The covariant components of the 
metric tensor 
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(12) 



together with the specified form of the energy momentum tensor T^^^ allows 
to derive the structure equations of a spherically symmetric star (Tolman 
Oppenheimer Volkoff equations) 



dP{r) 
dr 



(13) 



G 



(p(r)+p(r)) (-;:)+sy^ 



^^=Wp(r). (14) 
dr 

As the solution the mass-radius relation for white dwarfs is obtained. 
The analyzed model of the white dwarf has to include both the core of 
degenerate isothermal electrons and the atmosphere where the temperature 
changes with the radius. In this case the equation of hydrostatic equilibrium 
has to be supplemented by the transport equation 

(i(e^WT) _ 3 Kp L e^W ^^^^ 



dr 16a Airr"^ 



2Gm 



where e^^''^ is the redshift, L is the luminosity, k, is the opacity of the stellar 
material and a is the Boltzmann constant. After neglecting small corrections 



coming from the general relativistic form of the equation (15) and using 
the condition of isothermality T^, = const it is possible to determine the 
luminosity of a white dwarf. In the result the total radius of a white dwarf 
can be obtained from the relation 
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( Shapiro et al. 1983 ). In this relation i?* is the core radius and Jl is the 
molecular weight of the gas particles. The boundary conditions for the equa- 
tion ([l3|) are giving in the following form: P{r = 0) = Pc = P{pc) where 
Pc denotes the density at the center of the star which in this paper ranges 
from 10^ g cm~^ to 10^^ g cni~^. The condition P{r = R) = determines 
the star radius R. For such boundary conditions the value of the magnetic 
field and temperature are limited to the following values T = — 10^ K, 
7 = — 0.1, respectively. The results are presented in figures 1-3. 
Studying the properties of a star in the framework of the general-relativistic 
Thomas-Fermi model one can made the assumption that the temperature 
and chemical potential are metric dependent local quantities and the gravi- 
tational potential instead of Poisson's equation satisfies Einstain's field equa- 
tions. The energy momentum conservation T-j^'^ = which for spherically 
symmetric metric (??) can be written as 

dv 2 dP 

— = (17) 

dr P + p dr 

together with the Gibbs-Duhem relation and with the assumption that the 
heat flow and diffusion vanishes, (Israel 1976) give the condition 

^ = const. (18) 

This implies that the temperature and chemical potential become local met- 
ric functions 

T(r) = e~^«/2j.^ (;^9) 
;,(r) = e-^«/V. (20) 

where and /x* are constant. The first equation in (^) is the well known 
Tolman condition for thermal equilibrium in a gravitational field, (Tolman 
1934). Now, the Fermi-Dirac distribution function has to be redefined and 
is given in the following form 

1 ^ 1 

_ 1 

g(Vl + a:2 + 27n-z(r))/t(r)^_^ ' 

( Bilic et al. 1999). In the result instead of the OTV equilibrium equations 
together with the equation of state which has been calculated in the flat 
space-time one can derived three self-consistent equations (13,14,1^). These 



equations together with the local form of the equation of state are now func- 
tions of r. This is of particular importance in the case of strong gravitational 
fields. 
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3 Discussion 



In order to construct the mass-radius relation for white dwarfs the proper 
form of the equation of state have to be enumerated. For our needs we have 
chosen the plasma as the main constituent of the theoretical model of the 
white dwarf interior. All calculations have been performed at finite temper- 
ature and different from zero magnetic field and compared with those of zero 
temperature and without magnetic field. 

In Fig. 1 the applied forms of the equations of state for magnetic and non- 
magnetic white dwarfs have been presented. Dots indicate the location of 
critical points. The approximation of spherical symmetry can be used for 
densities grater than the critical one Pcrit- The inclusion of the magnetic field 
makes the equation of state stiffer. In the case of zero temperature model 
there are clearly visible Landau levels which became smear out when the fi- 
nite temperature is included. The figure in the subpanel shows the equation 
of state obtained for the low density range for the zero temperature case, Lan- 
dau levels are visible. In Fig. 2 the mass-radius relations for magnetic white 
dwarfs are presented. For each value of 7 the maximum mass has nearly 
the same value however the approximation of spherical symmetry which is 
connected with the existence of critical point determines the configuration 
with maximum radius. For the sake of completeness this figure includes also 
the mass-radius relation obtained for the equation of state of Hamada and 
Salpeter (Hamada et al. 1961). The solid line represents the nonmagnetic 
case which has been get with the use of the Hamada Salpeter equation of 
state. Dashed and dotted curves have been constructed for different values 
of the magnetic field strength. The differences between these curves are very 
small. The most important effect concerns the position of the critical points 
connected with the applied approximation of the spherical symmetry. In 
this approximation stable white dwarf configurations exist only for densities 
higher then pcrit- Thus the critical points on the curves determine configu- 
rations with maximal radius R. The stronger the magnetic field the smaller 
the white dwarf radius is. The most extended objects are obtained for the 
weakest magnetic fields. Using the Thomas-Fermi approximation one can 
obtain the result in which the star interior is no longer isothermal. The tem- 
perature profiles are presented in Fig. 3 for the fixed values of the central 
density pc- The temperature has been changed but not in a significant way. 
This is in agreement with the value of the gravitational potential which in 
the case of white dwarfs is much smaller than that obtained for neutron stars. 
The higher the central density the more visible are changes in temperature 
inside the star. For the central density = IQ^^gcm"^ the temperature dif- 
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[K] 


710*^ 


3.2398 10*^ 


310'^ 


2.211510'' 


3.5 lO^'^ 


2.2631 10*^ 



Table 1: The magnetic temperature Tb for the white dwarf with 7 = 0.01 
and temperature T = 0. 

ference in the center of the star equals 6 10^ K. The magnetic temperatures 
Tb for different values of the central density pc and different magnetic field 
strengths are presented in Tables 1, 2, 3. It has been found that Tb increases 
with the increasing value of the magnetic field and decreasing for increasing 
value of the central density. In Tables 4, 5, 6 the basic parameters of mag- 
netic and non-magnetic white dwarfs are collected. For the maximum mass 
configurations the atmosphere extent drop with the increasing strength of 
the magnetic field. The effective temperature rises with the increasing mag- 
netic field whereas the luminosity is altered insignificantly. All changes of 
white dwarf parameters are evident in the area of moderate densities. In this 
paper a magnetization A4 in the medium is not included. The inclusion of 
magnetization changes the perpendicular to the magnetic field component of 
a pressure and alters the value of the star radius. In the paper by Felipe it 
is stated that for positive magnetization the transversal pressure exerted by 
the charged particles in the magnetic field is smaller than the longitudinal 
one by the amount BAi^ (Felipe et al. 2002). Thus for n = (the ground 
Landau level) P± = 0. The vanishing of the pressure means that the insta- 
bility which leads to the gravitational collapse appear. However, considering 
the case of finite temperature it is necessary to take into account a fact that 
the temperature smearing out Landau levels causes the appearance of the 
residual pressure which supports tlie star against gravitational collapse. The 
instability criterion given by the condition 

n ~210-^^^ 

"-max — ^ '-^ 

sets in for densities equal 5.7 x 10^ g cm~^ for 7 = 0.01, 5 x 10^ g cm~^ 
for 7 = 0.05 and 1.5 x 10^ for 7 = 0.1. This instability leads to a cigar like 
object whereas in the case considered in this paper gives in the limiting case 
a toroidal star. 
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LOxio'n 




Figure 1: The equation of state for different cases of temperature and mag- 
netic field strength. Dots represent the position of critical points. The figure 
in the subpanel depicts the form of the equation of state constructed for the 
low density range with clearly visible Landau levels. 
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Figure 2: The mass-radius relations for magnetic white dwarfs obtained for 
different values of the magnetic field strength. The marked points represent 
the limiting values of radii for stable spherically symmetric configurations. 
The figure in the subpanel shows the mass-radius relations obtained for the 
radius range 1 — 1.25i?0/lOO. For comparison the nonmagnetic case is con- 
structed for the Hamada-Salpeter equation of state (Hamada et al. 1961). 
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Figure 3: The temperature profile for different central densities. 



Pc [g cm"^] 


Tb[K\ 


710*^ 


1.6199 10« 




1.1088 10« 


3.510^" 


1.131510''' 



Table 2: The magnetic temperature Tb for the white dwarf with 7 = 0.05 
and temperature T = 0. 



The numbers from 1 to 6 denote: 

1- the core radius of the white dwarf lOOi?^/i?0 

2- the total radius of the white dwarf lOO-R/i?© 

3- the extent of the atmosphere lOOAi?/i?0 

4- the effective temperature T^s [K], 

5- the critical density pcrit [9 cm^^], 

6- the strength of magnetic field 7 
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Pc [g cm ^] 


Tb[K] 


710*^ 


3.264710^ 


310^ 


2.2255 10« 


3.5 10^" 


2.2632 10^ 



Table 3: The magnetic temperature Tb for the white dwarf with 7 = 0.1 and 
temperature T = 0. 



1 


2 


3 


4 


5 


6 


1.083 


1.099 


0.0158 


13679 







1.089 


1.105 


0.0159 


15822 


4.610^ 


0.01 


1.087 


1.103 


0.0154 


15896 


3.6 10*^ 


0.05 












0.1 



Table 4: The parameters of the magnetic and non-magnetic white dwarf 
with the mass 0.74 Mq, central density pc = 710^g cm~^ and luminosity 
L = 0.0012 L„. 



1 


2 


3 


4 


5 


6 


0.807 


0.815 


0.00870 


24633 







0.810 


0.818 


0.00877 


30799 


4.6 lO-"" 


0.01 


0.809 


0.818 


0.00875 


30864 


3.6 10« 


0.05 


0.807 


0.816 


0.00872 


30991 


9.4 lO'* 


0.1 



Table 5: The parameters of the magnetic and non-magnetic white dwarf with 
the mass 1 Mq, central density Pc = 3 10^ g cmT^ and luminosity L = 0.0016 



1 


2 


3 


4 


5 


6 


0.130 


0.130 


0.00022 


940946 







0.117 


0.117 


0.00018 


1.6010*^ 


4.610^ 


0.01 


0.117 


0.117 


0.00018 


1.6110^ 


3.610*^ 


0.05 


0.111 


0.111 


O.OOOIG 


1.7810'' 


9.110'' 


0.1 



Table 6: The parameters of the magnetic and non-magnetic white dwarf 
with the mass 1.42 M©, central density pc = 3.5 10^°^ cm~^ and luminosity 
L = 0.0023 Lp,. 
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